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UNIMODULAR BINARY HIERARCHICAL MODELS 


DANIEL IRVING BERNSTEIN AND SETH SULLIVANT 

Abstract. Associated to each simplicial complex is a binary hierarchical model. We 
classify the simplicial complexes that yield unimodular binary hierarchical models. Our 
main theorem provides both a construction of all unimodular binary hierarchical models, 
together with a characterization in terms of excluded minors, where our definition of a 
minor allows the taking of links and induced complexes. A key tool in the proof is 
the lemma that the class of unimodular binary hierarchical models is closed under the 
Alexander duality operation on simplicial complexes. 


1. Introduction 

Associated to a simplicial complex C with ground set [m] and an integer vector d G 
is an integral matrix Ac,d- The hierarchical model associated to C, d is a log-linear model 
(i.e. discrete exponential family) whose design matrix is the matrix Ac,d- The Zariski 
closure of a hierarchical model is a (not necessarily normal) toric variety. As for all log- 
linear models, important relevant problems are to study properties of the lattice ker^ Ac,d, 
the polyhedral cone M>oAlc,d := {Ac,dX : a; > 0} and the semigroup NAc.d := {Ac,d^ '■ 
X > 0,a: integral }. For example, Markov bases of the lattice kery, A n d can be used for 
goodness-of-£t tests of the log-linear model associated to Ac,d (see |^). In the special 
case where d = 2 the vector of all twos, we call Ac,2 a binary hierarchical model. 

There are a number of “niceness” properties that an integral matrix A G could 
satisfy. Perhaps the strongest is unimodularity. 

Definition 1.1. The following are equivalent dehnitions for an integral matrix A to be 
unimodular. 

(1) The polyhedron P4= {x G : Ax = 6, x > 0} has all integral vertices for every 
beNA. 

(2) Suppose rankA = d. Then there exists a A such that the determinants of all d x d 
submatrices of A are 0 or ±A. 

(3) The entries in each circuit of A are either 0 or ±1. 

(4) The Graver basis of A consists of elements whose entries are either 0 or ±1. 

We elaborate on this dehnition in Section [2l A weaker property for the matrix A to 
satisfy is normality. 


Definition 1.2. A matrix A is normal if the semigroup NA is saturated, that is NA = 
ZA n M.>oA. 


In fact, every unimodular matrix A is normal. In previous work of Rauh and the 
second author , normality was identihed as a key property of hierarchical models to 
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be able to apply the toric fiber product construction to calculate a Markov basis. If 
Ac^d is unimodular, it is also easy to solve the integer programs that arise in sequential 
importance sampling j^. It is a major open problem to classify the normal hierarchical 
models. The special case where the underlying simplicial complex C is a graph, and where 
d = 2 was handled in [l^ where it was shown that a binary graph model is normal if and 
only if the graph is free of i ^4 minors. 

Our goal in the present paper is to classify the unimodular binary hierarchical models. 
In fact, this result is also an important first step in the classification of normal binary 
hierarchical models because a matrix A(74) of Lawrence-type is normal if and only if 
unimodular (this follows from [l^ Corollary 16 in light of the circuit definition of uni¬ 
modularity). So this provides a characterization of the normal binary hierarchical models 
that have a big facet (see Definition 13.141 and remarks that follow). 

Our main result is a complete structural characterization of the unimodular binary 
hierarchical models. On the one hand, we show that every unimodular binary hierarchical 
model can be built up from three basic families, together with three operations. As part 
of that proof we also give an excluded minor characterization of the unimodular binary 
hierarchical models. Here our notion of minor of a simplicial complex is one that is 
obtained by taking a sequence of vertex deletions or links of vertices. We caution the 
reader that our notion of simplicial complex minor does not directly generalize the notion 
of a graph minor, but it does generalize the notion of a matroid minor. Our excluded minor 
characterization involves an infinite number of minimally nonunimodular complexes. 

Our characterization of unimodular binary hierarchical models can be seen in analogy 
with the classic characterizations of regular matroids/totally unimodular matrices, which 
have both excluded minor characterizations and constructive characterizations. While 
there is clearly a close connection between our problem and the study of totally unimodu¬ 
lar matrices, our methods do not make use of Seymour’s [l^ nor Tutte’s [l^ classification 
of regular matroids. There are two reasons for this. The first is that the set of matroids 
underlying the design matrices of hierarchical models is not closed under the operations 
of direct sum, two sum, three sum, deletion, nor contraction. The second is that nor¬ 
mality, unlike (total) unimodularity, is not a matroidal property and we are interested in 
extending our methods to the normal classification problem [^. 

The outline of this paper is as follows. In the next section, we provide a detailed in¬ 
troduction to hierarchical models and the construction and interpretation of the matrix 
Ac,d- Section |3] details basic examples of unimodular complexes and constructions of new 
unimodular complexes from old ones. A main result here is Proposition 13.71 which shows 
that Alexander duality preserves unimodularity. Section H] gives the list of minimal non¬ 
unimodular binary hierarchical models. There are one infinite family and six sporadic 
minimal non-unimodular complexes. The section also provides a number of structural 
results on simplicial complexes that avoid this list of simplicial complexes. Section [S] 
provides the combinatorial description of the 1-skeleton of a simplicial complex that sup¬ 
ports a unimodular binary hierarchical model. That result is used as a stepping stone in 
a complex induction in Section |6] which gives the proof of the main theorem. Theorem 
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16.31 In Section [7] we explain how our characterization in the binary case can be used to 
make headway in the non-binary case. 


2. Preliminaries on Hierarchical Models and Unimodularity 

In this section we explain how to construct the matrix Ac,d associated to a hierarchical 
model. We then dehne some terms from Dehnition ll.il of a unimodular matrix and prove 
the equivalence between the four parts. 

Definition 2.1. Let d = {di ,..., dn) G Z"" be such that di>2 for each i. Then we dehne 
to be the vector space of all di x • • • x d„-way tables. That is, each u G is a table 
in n dimensions where the ith dimension has di levels. We denote entries in a table u 
by where (U, ■ ■ ■ ,in) ^ [di] x • • • x [d„]. We will use the shorthand i to denote 

Definition 2.2. Let C be a simplicial complex on ground set [n]. A facet of C is an 
inclusion-maximal subset F C [n] that is contained in C. We let facet (U) denote the 
collection of facets of C. 

Definition 2.3. Let d = (di,..., dn) G Z"" and F = {/i,..., /*,} C [n]. We dehne di? to 
be the restriction of d to the indices in F, i.e. di? = {df^,... ,df^). 

Definition 2.4. Let C be a simplicial complex on ground set [n] and let d = (di,..., d„) G 
Z"' be such that dj > 2 for each i. Then we have the linear map 

7rc,d : M'" ^ 0 

F£facet(C) 


dehned by 


/ 


■ ij e [dj]) !-)■ 0 


i^Gfacet(C) 




• j ^ [d/i] X • • • X [d/j,] 




F={fi,...Jk} \ iF=j 


The matrix .Ac,d denotes the matrix in the standard basis that represents the linear 
transformation vrc^d- 


We give an example of this construction. 


Example 2.5. Let C be the simplicial complex on ground set [3] = {1, 2, 3} with facet(C) = 
{{1},{2,3}}. The linear transformation 7rc,d maps a three-way tensor u = {uijk : i G 
[di], j G [d 2 ], k G [ds]) G to the direct sum of a one-way tensor and a two-way tensor: 




i,k 


G [di]) © C^^Uijk : i G [d 2 ], k G [ds]). 
i 
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Taking d = 2, the matrix Ac ,2 could be represented as follows: 


i 

0 

i 

0 

i 

0 

i 

0 

u 

1 

u 

1 

u 

1 

u 

1 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

VO 

0 

0 

1 

0 

0 

0 

V 


the horizontal line dividing the image space into the parts corresponding to the two facets 
ofC. 


It can be useful to be precise about how to write down the matrix Ac,d- The columns 
of the matrix Ac,d are indexed by elements of the set nje[n]['^j] rows are indexed 

by pairs {F,e) where F G facet(C) and e G The entry in position with row 

index {F,e) and column indexed by i G njG[n]['^j] is 1 if e = {ij : j E F), otherwise it is 
equal to 0. 

The binary case studied in this paper concerns the situation where all table dimensions 
are of size 2, that is d = 2 = (2, 2,..., 2); the vector of all twos. In this case, we will 
usually abbreviate notation as stated in the following definition. 

Definition 2.6. Let C be a simplicial complex on [n]. We define Ac ■= Ac, 2 - We say 
that C is unimodular if Ac is unimodular as a matrix. 

We now elaborate on Dehnition 11.11 To start, we give two definitions. 

Definition 2.7. Let A be an integer matrix. A nonzero element u G ker^A is called a 
circuit if its nonzero entries are relatively prime and there is no other nonzero element 
V G kergA such that supp(n) C supp(M). Here, supp(n) denotes the set of positions of 
nonzero elements of v. 


Definition 2.8. Let A be an integer matrix. A nonzero element u G ker^A is called 
primitive if there is no nonzero v G kerg A, v ^ u such that UiVi > 0 for all i, and 
ItjI < \ui\ for all i. The set of all primitive vectors in kerg A is called the Graver basis of 

A. 


Now we explain the equivalence between the four items in Definition 11.11 See [12 


Theorem 19.2 for the equivalence of (1) and (2). We can see that condition (3) of Definition 
ll.ll is equivalent to condition (2) because the circuits can be computed via a determinantal 
formula using Cramer’s rule [ij, p. 35]. Given the equivalence of (2) and (3), (3) and (4) 
are equivalent by [l^. Propositions 4.11 and 8.11. 

Questions about unimodularity of a matrix A are invariant under certain changes to 
the matrix. The following proposition gives three that we make tacit use of throughout 
the rest of the paper. 


Proposition 2.9. Let A be a unimodular integer matrix. Then the following are also 
unimodular matrices. 

(1) A matrix obtained by reordering the columns of A 
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(2) A matrix obtained by multiplying a column of A by —1 

(3) Any matrix that has the same rowspace as A 

(4) Any matrix formed by a subset of columns of A. 

Proof. Unimodularity of ([T]) and ([3]) are clear when considering (3) or (4) in Definition 
11.11 Unimodularity of ([2]) is clear when considering (2) in Definition ll.il If A' is a matrix 
obtained by extracting a set of columns of A then the circuits of A' are obtained by 
taking those circuits of A whose nonzero elements correspond to the columns of A we are 
extracting, and so dH) holds. □ 

3. Constructions of Unimodular Complexes 

In this section we describe some operations on simplicial complexes that preserve uni¬ 
modularity. In particular we show that unimodularity is preserved when passing to in¬ 
duced subcomplexes, to the Alexander dual of a simplicial complex, and to the link of a 
face of a simplicial complex. We can also build new unimodular complexes from old ones 
by adding cone vertices, ghost vertices, or taking a Lawrence lifting. We also construct the 
basic examples of unimodular complexes. These tools together go in to our constructive 
description of unimodular complexes. 

Proposition 3.1. All induced sub-complexes of a unimodular complex are unimodular. 

Proof. Let C be the induced subcomplex of C obtained by restricting to some vertex set 
F C [n]. Let B be the matrix obtained by taking the columns of Ac corresponding to the 
elements i = {ii,... An) £ U ~ ^ 3 ^ N \ ^et B' be the 

matrix obtained from B by removing rows of all zeros. Then B' = Ac- So Ac and B 
have the same Graver basis, which is a subset of the Graver basis of Ac- So unimodularity 
of Ac implies unimodularity of Ac- □ 

Next we will show that unimodularity is preserved under the Alexander duality opera¬ 
tion. 

Definition 3.2. Let C be a simplicial complex on [n]. The Alexander dual of C, denoted 
C* is dehned as: 

C* = {SG [n] : H \ 5 ^ C}. 

Note that if C has d faces, then C* has 2"’ — d faces and that the facets of C* are the 
complements of the minimal non-faces of C. 

Since we are now restricting attention to the binary case, where di = 2 for each i E [n], 
entries of vectors in are indexed by elements of {1,2}"' and entries of vectors in 
0/efacet(c) are indexed by pairs (S',}) where S G facet(C) and j G {1,2}'^. However, 
the statement and proof of Proposition l3.6l become cleaner if we replace {1, 2}" and {1, 2}'^ 
by (0,1}" and {0, l}"^, respectively. So we adapt this convention for the remainder of this 
section. We now introduce some notation that we use in Propositions 13.41 and 13.61 

Definition 3.3. Let S C [n] and j G {0, and i G {0, l}N\s. Then Cij denotes the 
element v G such that = 1 ii kt = jt when t G S' and kt = it when t G [n] \ S', 

and Tk = 0 otherwise. 
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We are now ready to give an explicit description of Ac*- 

Proposition 3.4. Let C be a simplicial eomplex on [n]. Let M he a matrix with the 
following set of columns: 

Cij : S' is a minimal non-face of C, i G {0, 

je{o,i}s 

Then Ac* = . 

Proof. We index the columns of M by pairs (S', i) where S is a minimal non-face of C, 
and i G {0,xhe rows of M are indexed by elements of {0,1}"'. By construction, 
M(k, (S', i)) = 1 if and only if k |[„]\5 = i. Similarly, when F is a facet of C*, .Ac* ((-S', i), k) = 
1 if and only kj^? = i. The result follows since S is a minimal non-face of C if and only if 
[n] \ S' is a facet of C*. □ 

In order to relate the unimodularity of Ac and Ac*, we need two propositions. The 
hrst is a standard result from the theory of matroid duality, though we could not hnd a 
precise reference. 

Proposition 3.5. Let A G have rank r. Let B G have rank n — r such 

that = 0. Then there exists a non-zero scalar A G K* such that for all S C [n] with 
#S = r, 

det(As) = ±Adet(-B[„]\5). 

Proof. After multiplying by an element of GL^(]K) and permuting columns, we can assume 
A has the form A = (A M) where A is an r x r identity matrix. Then we can suppose 
B has the form B = (— In-r)- Then any det(A 5 ) is a minor of M with row indices 
[r] \ S' and column indices S' \ [r]. The determinant det(-B[n]\ 5 ) is the same minor of M, 
up to a sign. □ 

Proposition 3.6. LetC he a simplicial complex on [n]. Then the following set of vectors 
spans ker(Alc)-' 

F = < (—: S' is a minimal non-face of C, i G {0, 

[j6{0,l}S 

Furthermore, if we view the entries of K as the columns of a matrix M', then we can 
multiply some set of rows and columns of M' by —1 and get M, as defined in Proposition 

H 

Proof. First, note that for any x E K, Acx = 0, so K C. ker(Alc). If we view the entries of 
K as the columns of a matrix M', we are done if we show that rank(M') = dim(ker(Alc))- 
We proceed by proving the second statement of the lemma. The hrst statement will follow 
since Proposition 13.41 implies rank(M) = rank(Alc*) and we know rank(.Ac*) = 2” — ffC = 
dim(ker(.Ac)) (@, Thm 2.6]). 

As before, we index the columns of M' by pairs (S', i) where S is a minimal non-face 
of C, and we index the rows of M' by the binary n-tuples. Let be a column of M'. 
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For any k G {0, !}"■, the entry ^ is nonzero if and only if k |[„]\5 = i - i.e. for each 

colnmn, the k for each nonzero entry has hxed k |[„]\5 = i. Now in this case, recall that 
= +1 if ||j||i is even, and = —1 if ||j||i is odd. Therefore, when ||i||i is 

even, = +1 if l|k||i is even and = —1 when ||k||i is odd. And when ||i||i is 

odd, k ~ when ||k||i is odd and k = “1 when ||k||i is even. So for a hxed 

colnmn the sign of a nonzero entry k depends only on the parity of k. So if 

we mnltiply all the rows M'k snch that ||k||i is odd by —1, the entries in each colnmn will 
all have the same sign. Then we can mnltiply all the negative colnmns by —1, to arrive 
at the matrix from Proposition 13.41 □ 

Now we are ready to show that Alexander dnality preserves nnimodnlarity. 

Proposition 3.7. Let C be a simplicial complex on ground set [n]. Then C is unimodular 
if and only if C* is unimodular. 

Proof. Consider the set K from Proposition 13.61 as a matrix of colnmn vectors. Since K 
spans the kernel of Ac, Proposition 13.51 implies that is nnimodnlar if and only if Ac 
is nnimodnlar. Then we can mnltiply the appropriate rows and colnmns of iF by —1 to 
get M as in Proposition 13.41 Since these operations preserve the absolnte valnes of fnll 
ra nk determinants, iF^ is nnimodnlar if and only if is. Proposition 13.41 says that 
Ac* = M^. □ 

Taking Alexander dnals and indnced complexes gives rise to another nnimodnlarity 
preserving operation which we now dehne. 

Definition 3.8. Let S' G C be a face of C. Then the link of S' in C is the new simplicial 
complex 

links(C) = {F\S :F eC and S CF}. 

When S = {n}, we simply write link^,(C) := hnk{^}(C). 

Note that we can obtain link 5 (C) by repeatedly taking links with respect to vertices. 
That is if S' is a face of C and ifS > 2 and n G S', then 

links (C) = hnk^(links\{^}(C)). 

Proposition 3.9. If C is a simplicial complex and S is a face of C, then hnks(C) = 

{c*\sy. 

Proof. By dehnition, we have: 

(C-\Sr = {iiC l„]\S:(ln|\S)\ii^C-\S}. 

Then we have the following chain of eqnivalences on some R F [n] \ S': 

^ ln|\((M\S)\fi)6C 
and since [n] \ (([n] \ S') \ i?) = i? U S', we have i? U S' G C. 


□ 
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Corollary 3.10. Let C be a unimodular simplicial complex on ground set [n]. Then for 
any face S ofC, link 5 (C) is unimodular. 

Proof. Proposition 13.91 implies that link 5 (C) can be obtained via Alexander duality and 
passing to an induced subcomplex. Unimodularity then follows from Propositions 13.II and 

[ATI □ 


Now we turn to operations for taking a complex that is unimodular and constructing 
larger unimodular complexes. If C has a vertex v that lies in each facet of C, we say that 
T is a cone vertex of C. The following proposition tells us that unimodularity is invariant 
under adding or removing a cone vertex. 

Proposition 3.11. Let C be a simplicial complex on [n]. We define C on [n + 1] to be 
the simplicial complex with the following facets: 

{F U {n + 1} : F is a facet of C}. 

Then Ac is unimodular if and only if Ac is. 


Proof. We will index the columns of Ac/ by the binary n + 1 tuples such that those with 
the n + 1 coordinate equal to 1 come before those with n + 1 coordinate equal to 2. Then 
we get the following block form: 


Ac 


Ac 

0 

0 

Ac 


The Graver basis of Ac is therefore {(m,0), (0,m) : u G Gr^}. Hence the Graver basis 
of Ac consists of 0, ±1 elements if and only if the Graver basis of Ac consists of 0, ±1 
elements. □ 


By induction, adding or removing multiple cone vertices from a simplicial complex does 
not affect unimodularity. We introduce the following notation to denote this. 

Definition 3.12. Let C be a simplicial complex on vertex set [n]. Then we dehne cone^(C) 
to be the simplicial complex on [n] U {mi, ..., Up} with the following facets 

{F U {til,..., Up} : F is a facet of C}. 

When p = 0, we dehne cone^(C) = C. 

Definition 3.13. For any matrix A G we dehne the Lawrence lifting of A to be the 
matrix 

[A 0\ 

A(A) = 1° ^1 ^ M(2^+bx2t 

where 0 denotes the s x t matrix of all zeroes and 1 denotes atxt identity matrix. 

By Theorem 7.1 in [l^, A(A) is unimodular if and only if A is unimodular. This gives 
rise to another unimodularity-preserving operation on simplicial complexes. 
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Definition 3.14. Let C be a simplicial complex on [n\. We define the Lawrence lifting of 
C to be the simplicial complex A(C) on [n + 1] that has the following set of facets: 

{N} U {F U {n + 1} : F is a facet of C}. 

In this case, we refer to the facet [n] as a big facet. 

If a complex C on [n + 1] has big facet [n], then C = A(link„_|_i(C)). In particnlar, 
C = A('D) for some complex F if and only if C has a big facet. 


Proposition 3.15. The simplicial complex C is unimodular if and only if K{C) is uni- 
modular. 


Proof. Note that by construction A(^c) = •Aa(c)- Hence Theorem 7.1 in 
proposition. 



implies the 

□ 


If a simplicial complex C can be expressed as A(C'), then we say that C is of Lawrence 
type. A simplicial complex on n vertices is of Lawrence type if and only if it has a facet 
containing n — 1 vertices. We will refer to any facet of C that has n — 1 elements as a big 
facet. 

We now give our hnal unimodularity preserving operation. 


Definition 3.16. Let C be a simplicial complex on ground set [n]. Let G^{C) denote the 
same simplicial complex but on ground set [n-\-p]. Note that the vertices n + 1,..., n + p 
are not contained in any face of G^{C). In this case we say that n + 1,..., n + p are ghost 
vertices. When p = 1 we drop the superscript; i.e. we just write G{C). 

Proposition 3.17. A simplicial complex C is unimodular if and only ifG^{C) is unimod¬ 
ular. 


Proof. This is true because 

aIgp(c) = (A: • • • A) • n 

Note that adding a ghost vertex to a complex is Alexander dual to taking the Lawrence 
lifting of a simplicial complex, that is G{C)* = A(C*). 

We now state a useful fact about the interaction of these operations. 

Proposition 3.18. The operation cone^(-) commutes with the operations of taking links, 
duals, induced sub-complexes and adding ghost vertices. 

For any vertex v of C, we let C \ n denote the induced subcomplex on C \ {n}. Then for 
vertices v ^ u oi C, the operations ■ \ v and link„(-) commute. So if F can be obtained 
from C by applying a series of deletions, and taking links, then we can write 

F = linkK(C \ S) 

where F E C and S, R are a subset of vertices of C such that S' fl F = 0. This gives rise 
to the following dehnition: 

Definition 3.19. Let C be a simplicial complex. Then F is a minor of C if F = hnk^(C\S') 
where and S', R are subsets of vertices of C such that S fl F = 0, and F is a face of C. 
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We warn the reader that although it is natural to view simplicial complexes as a gen¬ 
eralization of graphs, our dehnition of simplicial complex minor does not generalize the 
usual notion of graph minor. Alternatively, a simplicial complex can be seen as a general¬ 
ization of a matroid, and our dehnition of minor is a generalization of the usual notion of 
a matroid minor [^. This dehnition of simplicial complex minor is useful for our purposes 
because of Proposition I3.201 

Proposition 3.20. If C is a unimodular simplicial complex, then every minor of C is 
unimodular. 

Proof. This follows immediately from Propositions 13.11 and 13.101 □ 

The fundamental example of a unimodular simplicial complex is the disjoint union of 
two simplices. Here a simplex A„ is the simplicial complex on an n -f 1 element set with a 
single facet consisting of all the elements. We denote the irrelevant complex {0} by A_i 
and the void complex {} by A_ 2 . 

Proposition 3.21. LetC = AmLJA„ he the disjoint union of two simplices, form,n > 0. 
Then C is unimodular. 

Proof. The matrix Ac in this case is the vertex-edge incidence matrix of a complete 
bipartite graph with 2"^^^ and 2”+^ vertices in the two parts of the partition. Such 
vertex-edge incidence matrices are examples of network matrices and are hence totally 
unimodular [H, Ch. 19]. □ 

A second family of fundamental examples comes from taking the duals of the disjoint 
union of two simplices. We use Ilm,n to denote the dual of a disjoint union of an m-simplex 
and an n-simplex, i.e. 

D •= lA U A P 

We close this section by giving a workable description of Dm,n- We can divide the vertices 
of Dm,n into disjoint sets M, N such that M contains the vertices of the A^ in „ and 
N contains the vertices of the A„ in Then, the facets of Dm,n are precisely the 

subsets of M U A that leave out exactly one element of M and one element of N. Notice 
that the complexes induced on M, N are d/S.rn and c?A„ respectively. Also, note that for 
any v e M, \mk^{Dra,n) = Dm-pn and for any v e N, hnk„(T)m,n) = Dm,n-i- 

4. /9-avoiding Simplicial Complexes 

Part of our main result is a forbidden minor classihcation of unimodular simplicial 
complexes. In this section, we identify these forbidden minors and prove various properties 
about the complexes that avoid them. 

Proposition 4.1. The following simplicial complexes are minimal nonunimodular sim¬ 
plicial complexes: 

(1) P 4 , the path on 4 vertices 

( 2 ) Oq, the boundary of the octahedron or its dual Og 

(3) Ji, the complex on {1,2, 3,4, 5} with facets 12,15,234,345 or its dual J*, with 
facets 134,235,245. 
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(4) J 2 , complex on {1, 2, 3,4, 5} with facets 12, 235, 34,145 

(5) For n > 1, dAn U {v}, the disjoint union of the boundary of an n-simplex and a 
single vertex. 



Figure 1 . The complexes P4, Ji, J^, and J2. 


Note that P 4 , J 2 and dAn U {n} are isomorphic to their own duals. 

Proof. We can check that the complexes P 4 , Oq, Oq, Ji, J*, and J 2 are not unimodular by 
using the software 4ti2 [l| to compute the Graver basis and looking for entries that are not 
0, ±1 In the case of Og and Og*, these are too large to compute the entire Graver basis. 
However, selecting sufficiently large random subsets of the columns produced Graver basis 
elements of the desired form. Gomputations are available on our website [^. For the 
infinite farnily dAn U {n} where n > 1, examples of a non-squarefree Graver basis element 


appear in [l5j. These results show that these examples are not unimodular. To see that 


they are minimal, note that every subcomplex obtained by deleting a single vertex or 
taking the link at a vertex produced a unimodular complex in all cases. □ 


Definition 4.2. A simplicial complex C is fd-avoiding if it does not contain any of the 
complexes from Proposition 14.11 as minors. 

Since none of the complexes from Proposition 14.11 is unimodular. Proposition 13.201 im¬ 
plies that if C is unimodular, then C is /9-avoiding. The converse of this is our forbidden 
minor classification of unimodular complexes which we prove in Theorem 16.31 

We now give some technical results about /9-avoiding complexes. Before beginning, we 
remind the reader that cIAi U{n} is an independent set on 3 vertices, and so no /9-avoiding 
complex can have an independent set of size 3. 


Proposition 4.3. Let C be a fd-avoiding simplicial complex. Then C* is also (d-avoiding. 

Proof. The list of prohibited minors of /9-avoiding complexes is closed under taking duals. 
The proposition follows when we note that if P is a minor of C, then P* is a minor of C*. 
This is true because if V* = link 5 '(C* \ R), then two applications of Proposition 13.91 give 

V* = {{c* \ R)* \ sy 

= (link^(C) \ 

and therefore P = link/^(C) \ S'. So "D is a minor of C. □ 

Proposition 4.4. Let C be a (d-avoiding simplicial complex that has G 4 induced. Then 
the complex induced on the non-ghost vertices of C is cone^^Cy for some p. 
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Proof. Assume C has no ghost vertices. Let ui,...,U 4 be the vertices that induce C 4 
in C. If these are the only vertices, we are done, so let v be another vertex in C. Let 
P denote the complex induced on v,ui,... , 114 . Any minor of P is a minor of C, so V 
must also be /9-avoiding. Vertex v cannot be disconnected from ui,... ,U 4 since otherwise 
P has an independent set of size 3. Furthermore, v must connect to Ui,... ,U 4 ^ with a 
triangle. Otherwise if v connected to Ut with an edge, then Ui has edge degree 3 and so 
\mk^{T>) is an independent set on 3 vertices. There are 4 possible triangles that v can 
join with wi,..., M 4 in "D so that ui,..., M 4 induce O 4 . If v is only in one such triangle, 
wlog {n, Ml, 112 }, then {m,mi,M 3 ,M 4 } is an induced P 4 . The two non-isomorphic complexes 
on 5 vertices that have two triangles and an induced C 4 , are Ji and J 2 . If P has 3 of the 
possible triangles, then link^(P) is P 4 . So V must have all four - i.e. it must be cone^( 04 ). 

Any vertices ^ {mi, M 2 , 114 } must connect to each other. Otherwise, the induced 

complex on v,v',Ui, ... ,M 4 is Og which is not unimodular. Now we denote the vertices 
of C that are not any of the m^s as Vi,...,Vk. Since Vk is a cone point over the 
on mi,...,M 4 and since Vk connects to all Vj,j < fc, hnk^,^(C) is a /9-avoiding complex 
on Ml,..., Vk-i, Ml,..., M 4 that has C 4 induced on mi, ..., M 4 . So, by induction on the 
number of vertices, link^j,(C) is an iterated cone over O 4 , cone^“^( 04 ). But also the 
induced complex on mi, ..., Vk-i, Mi, ..., M 4 is cone^“^(C' 4 ) by induction. This implies that 
C = cone^ {C 4 ). □ 

Proposition 4.5. LetC be a (^-avoiding complex that has Dm,n induced for some m,n > 1. 
Then the complex induced on the non-ghost vertices of C is cone^(Pm,n) for some p. 

Proof. Assume C has no ghost vertices. Just as in the remarks at the end of section [3l 
we divide the vertices of Dm,n into disjoint sets M, N such that M contains the vertices 
of the Am in N contains the vertices of the in Now, we proceed by 

induction on m and n. Note that Pi^ = so the base case is handled by Proposition 
14.41 Assume that, without loss of generality, m > 2 . 

Let V E C \ Dm,n- The vertex v must connect to some m G M to avoid inducing 
dAm U {m}. Then, link„(C) contains v and Pm-i,n, and so by induction, the complex 
induced on v and the vertices of Pm-i,n in hnk„(C) is a cone over Pm-i,n- This means 
that V is in every facet that contains m. Since m > 2, m is connected to every other vertex 
in Dm,n and thus v is attached to every vertex in Pm,n- Now we apply the same argument 
to each of the vertices in the set M to see that v is in every facet that contains any vertex 
in set M. Every facet of Pm,n contains some element of M. So this implies that the 
induced complex on v and Dm,n must be the cone over Dm,n- 

Now assume m, v' are both vertices in C \ Dm,n- If they were not connected by an edge, 
then the induced complex on m, v' and Dm,n contains a minor which is isomorphic to 
Oq. This means that v, v' are connected by an edge. Taking the link hnk^(C) produces a 
smaller complex with an induced Pm,n hence it must be a cone by induction on the number 
of vertices not in the Pm,n- This reduces us the case where hnkt,(C) is cone^“^(Pm,n) and 
C \ M is cone^“^(Pm,n) which implies that C is cone^(Pm„n)- D 
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Proposition 4.6. Let C be a (3-avoiding simplicial complex on vertices ui,... ,Uk+i,v 
such that the complex induced on {mi, ... ,Mfc+i} is dAk- Then v must be in a k-simplex 
with some subset of the Ui’s. 

Proof. We proceed by induction on k. For the base case k = 1, note that c?Ai is two 
isolated points. In this case, v must connect to Ui or U 2 as a 1-simplex (an edge) to avoid 
inducing dAi U {n}. 

Now assume k > 1. The vertex v must attach to some Ui to avoid inducing dAk U 
{n}. Then, hnk(Mj) has v and dAk-i, so by induction, v must form a k — 1 simplex 
with some collection ui,... ,Ui,... ,Uj,..., Uk+i- So in C, v must be in a fc simplex with 
Ml, . . . , Ujj . . . , Mfc-i-i. IT 

Proposition 4.7. Let C be a (3-avoiding simplicial complex on m + n + 2 vertices that 
does not have a facet of dimension m + n. Assume Dm,n ^ C. Then C = Dm,n- 

Proof. Note that C* C U A„ because Dm,n ^ C. Since C does not have a facet of 
dimension m-\-n, C* has no ghost vertices. Since C* has no ghost vertices, if C* C AmU A„, 
then C* would have an induced {m} U dAk for some k >1. In this case C* would not be 
/9-avoiding and so by Proposition 14.31 neither would C. □ 

5. The 1-Skeleton of a /9-avoiding Complex 

In this section we prove Lemma 15.31 which gives a complete characterization of the 1- 
skeleton of a /3-avoiding simplicial complex. This is a crucial technical lemma in the proof 
of Theorem 16.31 We start with a technical proposition about graphs. 

Proposition 5.1. Let H be a connected graph that avoids and P 4 as induced subgraphs. 
Then H is a complete bipartite graph. 

Proof. Let u G V{H), let N{u) denote the neighbors of u, and let M{u) denote the non¬ 
neighbors of u (this set includes u). The bipartition of the vertices of LI will be M{u) and 
N{u). Let V G M{u)\{u}. Since H is connected, there exists a path u = ui,U 2 ,... ,Uk = v. 
Assume k is minimal. Since v ^ u, k > 1. We cannot have k = 2 since u,v are non¬ 
neighbors. We cannot have fc = 4, since in order to avoid an induced P 4 , we would need an 
edge (Mi,Mi+ 2 ) inducing a P’s, or an edge (m,m) contradicting that m,m are non-neighbors. 
If A; > 5, there must exist an edge (mi, M 4 ) to avoid an induced P 4 , contradicting minimality 
of k. So we have k = 3. So for any v G M{u) \ {«}, there exists a path u, a, v. 

Now we show that H is bipartite with bipartition M{u) and N{u). It is clear that 
N{u) is an independent set of vertices, for if m,m; G N{u) had an edge between them, 
there would be a P 3 on u, v, w. Now we show that M{u) is an independent set of vertices. 
Assume w,v & M{u). If either w,v is u, there is no edge between them by dehnition of 
M{u), so assume w,v ^ u. Then by the above, we have paths u,a,v and u,h,w. So 
an edge (m, w) would induce the 5-cycle u, a, v, w, b, u. This 5-cycle must have a chord to 
avoid inducing P 4 , but any chord in a 5-cycle induces a K^. 

Now we show that H is complete bipartite. Let x G M{u) and y G N{u). Since H is 
connected, there is a path x = Wi, ...,ma: = U- We may without loss of generality assume 
k <3 since otherwise we could shorten the path using the edge required to avoid 
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a P 4 . The sets M{u) and N{u) are disjoint, so fc 7 ^ 1. If fc = 3, then there is an induced 
P 4 X, U 2 , y, u. So k = 2, and so (x, y) is an edge. □ 

For a graph G we dehne the complement graph G‘^ on the same set of vertices such that 
(u, v) is an edge of G^ if and only if (n, v) is not an edge in G. Now we can use Proposition 
15.11 to give a strong restriction on the structure of G'^ whenever G is the 1 -skeleton of a 
/3-avoiding complex. 

Proposition 5.2. If G is the 1-skeleton of a fl-avoiding simplicial complex C, then each 
connected component of G^ is complete bipartite. 

Proof. Note that {n} UclAi consists of three disconnected vertices, its complement graph 
is K^. The path P 4 is its own complement. Since C is /3-avoiding, G avoids P 4 and an 
independent set of size three as induced subgraphs. So each connected component of 
G'^ avoids P4 and as induced subgraphs. Proposition 15.11 therefore implies that each 
connected component of G‘^ is complete bipartite. □ 

Now we are ready to characterize the 1-skeleton of a /3-avoiding simplicial complex. 

Lemma 5.3. Let C be a ft-avoiding simplicial complex and let G denote its 1-skeleton. 
Then G is one of the following 

(a) Kn 

(b) Two complete graphs glued along a (possibly empty) common cligue 

(c) The iterated cone over a i-cycle. 

Proof. By Proposition 15.21 we know that each connected component of is complete 
bipartite. Let H denote the induced subgraph of G'^ that removes all isolated vertices. If 
H is empty, then G^ is an independent set of vertices and therefore G = Kj^. So assume 
H is nonempty. We claim that if H is neither Km,n nor K 2 U K 2 for m,?7, > 1, then H 
induces either K 2 U K 2 U K 2 or P 3 U K 2 . To prove the claim, hrst assume that H avoids 
K 2 U K 2 U K 2 . Since H is not Km,n and has no isolated vertices, this implies that H has 
exactly two components. Since H is not K 2 U K 2 , some connected component of H has 
at least three vertices. So we have a P 3 induced by this component, and K 2 induced by 
the other component and so the claim is proven. 

The complement graphs of K 2 U K 2 U K 2 or P 3 U K 2 are shown in Figure [2J Both 
graphs have G4 induced, but neither is the 1 -skeleton for an iterated cone over G4 since 
neither graph has a suspension vertex. Proposition 14.41 therefore implies that neither is 
the 1 -skeleton of a /3-avoiding simplicial complex and so G may not induce either. So G^ 
may not induce K 2 U K 2 U K 2 nor P 3 U K 2 and therefore neither may H. The claim then 
implies that H must be either or K 2 U 7^2 • Assume G^ has p isolated vertices. If 
H = Km,n then G is a K^+p and a Kn+p glued along a common Kp. li H = K2 U K2 then 
G is an iterated cone over a 4-cycle. □ 

As we see in the following proposition, the 1 -skeleton of C completely determines C 
when it is obtained by gluing two complete graphs along an empty clique. 

Proposition 5.4. Let C be a fl-avoiding simplicial complex such that the 1-skeleton of C 
is the disjoint union of cligues Km and K^. Then C = Am U A„. 








UNIMODULAR BINARY HIERARCHICAL MODELS 


15 




Figure 2. The complement graphs of P 3 U K 2 and P 2 LI iF 2 LI K 2 

Proof. Let Vi,..., be vertices of the K^- If {ui,..., Vk\ is a minimal non-face of C, then 
if we let M be a vertex in the then the complex indued on ui,..., u is dAkU{u}. □ 

6 . The Main Theorem 

The goal of this section is to give a proof of Theorem 16.31 which gives a complete 
characterization of the unimodular binary hierarchical models. 

We begin this section by defining a nuclear complex. A nuclear complex is a complex 
that can be obtained from a disjoint union of two simplices by adding cone vertices, adding 
ghost vertices, taking Lawrence liftings and taking Alexander duals. Since L A„ 
is unimodular and these operations all preserve unimodularity, nuclear complexes are 
unimodular. Part of our main result is the converse - unimodular complexes are nuclear. 

Definition 6.1. A simplicial complex C is nuclear if one of the following is true 

(1) C = A(P) where P is nuclear 

(2) C = G(V) where V is nuclear 

(3) C = cone^(Am L A„) for p,m,n >0 

(4) C = cone^(Dm,n) for p > 0 and m, n > 1 

(5) C = Afc for k > — 2 . 

Every nuclear complex C can be constructed by applying the operations cone^(-), 
and A(-) to a complex T) where P is of the form A^ L A„, Dm,ni or A^. We refer to P as 
the nucleus of C. 

Note that Dm,o has a ghost vertex for all m. This is why we have m, u > 1 in (|1]). We 
note that the collection of nuclear complexes is closed under Alexander duality. 

Proposition 6.2. If C is nuclear then so is C*. 

Proof. Possibilities ([T]) and ([ 2 ]) are dual to each other, as are ([3]) and (jl]). On k vertices, 
Afc and A _2 are dual to each other. □ 

We now state our main result. 

Theorem 6.3. Let C be a simplicial complex. Then the following are equivalent 

(1) C is unimodular 

(2) C is (3-avoiding 

(3) C is nuclear. 
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We defer the proof of Theorem 16.31 until the end of the section, but we give a roadmap 
here. It is immediate from results in previous sections that nuclear complexes are uni- 
modular and that unimodular complexes are /3-avoiding. 

We show that any /3-avoiding complex C is nuclear by choosing a particular vertex v 
of C and using induction on the number of vertices to conclude that link^(C) is nuclear. 
From here, we have hve cases to consider - one for each of the ways that linkj,(C) can 
be nuclear. In four of these cases, it is relatively easy to show that C is nuclear. The 
difficulty lies in the case where hnkt,(C) = Dm,,n, which is handled in Lemma 16.81 The 
proof of Lemma 16.81 is split into two main parts and the second part is further split into 
six cases. Proposition 16.41 helps with the hrst part. Propositions 16.5116.61 and 16.71 handle 
the hard cases in the second part. 

Proposition 6.4. Let C be a [^-avoiding simplicial complex. Assume that all proper 
minors of C are nuclear. Assume C has a vertex v such that C\v = cone^(Am U A„) 
with m,n > 1 and p > 1. Assume C is C^-free. Then C = cone^’''^(Am U An), C = 
cone^’(Am+i U A„), or C = coneP(Am U A„+i). 

Proof. We induct on p. For the base case, take p = 1. If we let u denote the cone vertex 
in C \ n then v must connect to one of the simplices in the Am U A„ G C\u. Otherwise 
we have an independent set with three vertices induced in C \ n and this complex is not 
nuclear. If v connects to only one such simplex, then v must also connect to u to avoid 
inducing P 4 which is not nuclear. If v connects to both such simplices, and v does not 
connect to u, then C has O 4 induced which contradicts our hypothesis. So assume v 
connects to u and to at least one simplex. Since n is a cone vertex in C \ n, link^t(C) has 
no ghost vertices and by hypothesis, it is nuclear. Furthermore, link„(C \v) = Am U A„ 
and so hnk^(C) is either A^+i U A„ or cone^(Am U A„) (note that these are the only 
nuclear complexes that become disconnected upon removing a vertex). In the hrst case, 
C = cone^(Am+i U A„) and in the second case C = cone^(Am U A„). 

Now assume p > 2. Let ui,... ,Up denote the cone vertices of C \ n. Let M, N denote 
the vertex sets of the A^, A„ respectively. Then v must connect to at least one of M and 
N to avoid inducing an independent set on three vertices. Furthermore, v must connect 
to each Ui - ii v only connects to one of M or N, then this is required to avoid P 4 , and 
if V connects to both then we need this to ensure that C is C' 4 -free. Furthermore, the set 
{«!,... ,Up,v} is a facet of C. This is clear when p = 1 since otherwise v doesn’t connect 
to Ml. Induction on p gives that {ui,... ,Mp_i,n} is a facet of linkup(C) and therefore 
{ui,..., Up, v} is a. facet of C. 

Now, we can see that linkup (C) has no ghost vertices and is C 4 -free. Furthermore, 
linktip(C) \ V and C \ {up, v} are both equal to cone^“^(Am U A„). So by induction on p, 
each of linkup (C) and (C \ Up) can either be cone^(Am U A„) or cone^“^(Am+i U A„). If 
linktip(C) = C\up, then Up is a cone vertex in C. So if both equal cone^(Am U A„) then 
C = cone^’''^(Am U A„). If both equal cone^“^(Am+i U A„) then C = cone^(Am+i U A„). 

Since link^jp(C) C C\up, it is impossible to have linkup(C) = cone^(Am U A„) and 
C \ up = cone^“^(Am+i U A„). It is also impossible to have C \ up = cone^(Am U A„) 
and linkup (C) = cone^“^(Am+i U A„). For if this were the case, then v would be a 
cone vertex in C \ Wp and it would be part of the A^+i in linkup (C). Since n > 1, 
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7 ^TV > 2, SO given bi,b 2 E N and a E M the complex induced on a, 6 i, 62 , Up, v has facets 
{n, a, Up}, {n, 61 , 62 }, {up, &i, & 2 }- This is J}, and so C is not /T-avoiding. □ 

Proposition 6.5. Let C he the simplicial complex on q vertices with a vertex v such that 

linkt,(C) = Dm,n with m,n > 1 and (C \ v)* = (A^' U A„) with 

and 0 < n' < n. Unless m' = —1 and n' = n, C has a proper minor that is not nuclear. 

Proof. Let M, N denote the sets of vertices such that in (hnkt,(C))*, the complex induced 
on M is Am and the complex induced on N is A„. Then some vertex a; G A is a non-ghost 
in (C \ n)*. We show that link 2 ;(C) is not nuclear. 

First we describe the facets of linkj.(C). We claim that the facets of hnkj;(C) are precisely 
the sets containing A — 1 — 2 vertices, where the omitted vertex pair is one of the following 

(1) some a E M and some b E N \ {x} or 

( 2 ) V and some vertex g that is a ghost in (C \ n)*. 

We now prove the claim. The facets of hnka,(C) that contain v are the sets whose com¬ 
plements are pairs of the first type since linki,(C) = Dm,n- The facets of hnkj,(C) that do 
not contain v are the sets whose complements are minimal non-faces of (C \ v)* that do 
not become faces of linka;(C) upon adding v. These minimal non-faces of (C \ n)* are the 
ghost vertices. Therefore the facets of linka;(C) that do not contain v are the sets of the 
second type. 

We now show that linka;(C) is not nuclear. The complex hnka,(C) cannot be A(D) for 
any complex V since no facet lacks fewer than two vertices. Nor can hnka,(C) have any 
ghost vertices - m > 1 implies ffM > 2 and so every vertex in M is in a facet of the hrst 
type, every vertex of A is in a facet of the second type, and v is in every facet of the hrst 
type. It is clear that no vertex is in every facet, so hnkj,(C) cannot be cone^(Afc;) nor 
cone^(Afc U A;) for any p > 1. It is clear that hnka,(C) has more than 2 facets, so hnka,(C) 
cannot be A^ U A;, nor A^. 

The only case left to check is hnkj,(C) = D^ i. If this is the case then we can partition 
the vertices of hnka,(C) into disjoint sets Ci^C^ such that the facets of hnka,(C) are pre¬ 
cisely the subsets of Ci U C 2 that omit exactly one from each Ci and 6 * 2 . 

Assume n' = n and m' > 0. Let g E M he a. ghost vertex and without loss of generality 
assume g E Ci. Facets of the hrst type imply M C Ci. Facets of the second type imply 
V E C 2 . Since m' > 0, there exists some a E M that is not a ghost in {C\v)*. Snd so 
(Cl U C 2 ) \ {a, v} is a facet of \mkx{C). But this does not ht our earlier description of the 
facets of linkj;(C). 

Now assume — 1 < m' < m and 0 < n' < n. Then there exist g' E M and g” E A that 
are ghosts in (C \ n)*. But this contradicts our earlier description of the facets of hnka,(C) 
- facets of the hrst kind imply that g', g" are in diherent CjS and facets of the second kind 
imply that they are in the same Cj. □ 

Proposition 6.6. LetC be a simplicial complex on q vertices. Let v be a vertex in C such 
that linkt,(C) = Dm,n with m,n> 1. Assume {C\v)* = G'^~^^{An). Then C = Dm,n+i- 

Proof. Let M, A denote the sets of vertices such that in (hnkj,(C))*, the complex induced 
on M is Am and the complex induced on A is A„. Then all the vertices of M and none 
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of the vertices of N are ghosts in (C \ n)*. We claim that each facet of C contains N — 2 
vertices, and the omitted pair is one of the following 

(1) some a & M and some b & N 

(2) V and some a G M. 

To see this, note that if a facet contains v, then it mnst be of the hrst form since 
link„(C) = Dm,n- If a facet does not contain n, then its complement mnst be n, along with 
a minimal non-face of (C \ n)*. These minimal non faces are precisely the vertices in M. 
So C = Dm,n+i where M and iV U {n} are the sets of vertices for the and A^+i in C* 
respectively. □ 

Proposition 6.7. Let C be a simplicial complex on q vertices. Let v be a vertex in C 
such that link.u(C) = T*m,n with m,n > 1 and has no ghost vertices. Let M,N denote the 
sets of vertices such that in (link„(C))*, the complex induced on M is Am and the complex 
induced on N is A„. Assume {C\v)* = where V is a nuclear complex on vertex 

set N with more than one facet. Then C has a proper minor that is not nuclear. 

Proof. Let x G M. We show that hnk 3 ;(C) is not nnclear. We start by showing it is not 
Lawrence type. Let F be a facet of hnk 3 ,(C) that does not contain v. Then F U {a:} is 
a facet of C \ n so it must also lack some a & M \ {x}. If F is a facet of link 3 ;(C) that 
contains v, then F U {x} \ {n} is a facet of hnk^,(C) and therefore lacks some a E M and 
b E N. So linka,(C) is not A(F) for any complex F. 

We can see that linka,(C) has more than two facets - since m,n > 1, there are at least 
two that include v and since {C \ v)* = there is at least one that does not. 

link 2 ,(C) cannot be cone^(Afc U A/) for any p > 0. 

All of the vertices of M are ghost vertices in (C\n)*, and (C\n)* has a minimal non-face 
S F N with at least two vertices. This implies that C contains a facet F such that M <E F, 
there exist 61,62 £ Af such that 61,62 ^ F and v ^ F (since linkt,(C) = Dm,n)- Facets in 
link 3 .(C) that contain v have A — 1 — 2 vertices - they lack some a E M \ {x} and some 
b E N. Furthermore, F \ {x} is a facet of linka;(C) containing at most A — 1 — 3 vertices. 
This means that hnk 3 ;(C) is not pure, and therefore not cone^{Dk^i) for any p > 0. So 
link 2 ;(C) is a proper minor of C that is not nuclear. □ 

Lemma 6.8. Let C be a /3-avoiding simplicial complex on q vertices. Assume all proper 
minors of C are nuclear. If there exists a vertex v of C such that linkt,(C) = cone^(Fm,n) 
for p > 0 and m,n >1, then C is nuclear. 

Proof. We split this into two cases. For the hrst case, assume p > 1. Here, Propositions 
13.91 and 13.181 give C* \ v = cone^(Am U A„). By Proposition 14.31 C* is /^-avoiding so in 
C*, V must connect to one of the simplices in the A^ U A„ C C* to avoid inducing an 
independent set on three vertices. Now let m be a cone vertex. If v connects to only 
one such simplex, then v must also connect to u to avoid inducing F 4 . If n connects to 
both such simplices, and v does not connect to u, then C* has F 4 induced. In this case 
Proposition 14.41 implies C* = cone^(F 4 ) and so C = cone*^(Ai U Ai). So we can assume 
V connects to at least one simplex, and to all cone vertices in C* \ n. From this we can 
see that C* has no induced G 4 . Since every proper minor of C* is nuclear. Proposition 16.41 








UNIMODULAR BINARY HIERARCHICAL MODELS 


19 


implies C* = cone^+^(Am U A„) or C* = cone^(Am+i U A„). Proposition 16.21 then implies 
that C is nuclear. 

For the second case, assume p = 0. So linkt,(C) = Dm,n with m,n> 1. Then Dm,n ^ 
C\v and so {C\v)* C A^ U A„. Let M, N denote the sets of vertices of C such that 
in (link^(C))*, the vertices of the A^ are M, and the vertices of the A„ are N. The 
non-ghost vertices of (C \ n)* must form a nuclear complex, and the only disconnected 
nuclear complexes are of the form A^ U A*. From this it follows that {C \ v)* must be one 
of the following forms (without loss of generality) 

(1) A^U A„ 

(2) U A„/) with —l<m'<m and 0 < n' < n, but not both m' = — 1 
and n' = n 

(3) 

(4) where D is a nuclear complex on N with more than one facet 

(5) G™+^+2('|0}) 

( 6 ) 

We handle each possibility separately. 

Case 1. Assume (C \ v)* = A^ U A„. Then (C \ v)* has no ghost vertices, and so C \ n 
is not of Lawrence type. So, by Lemma ITTl C\v = Dm,n and so C = cone^{Dm,n) with v 
as the cone point. 

Case 2. Assume (C \v)* = {^m' LI A„/) with — 1 < m' < m and 0 < n' < n 

but not both m' = —1 and n' = n. By Proposition 16.51 C has a proper minor M that is 
not nuclear. By induction, M is not /^-avoiding and so neither is C. 

Case 3. Assume {C\v)* = G'^^^{An). Then Proposition 16.61 implies that C = Dm,n+i- 

Case 4. Assume {C \ n)* = G'^^^{D) where D is a nuclear complex on N with more 
than one facet. By Proposition 16.71 C has a proper minor M that is not nuclear. By 
induction, M is not /^-avoiding and so neither is C. 

Case 5. Assume {C \ v)* = G'^^'^^‘^{{^}). Then C\v = d^m+n+i- But link^C = Dm,n 
and so v is not in any m -|- n — 1 simplex with the vertices of M and N. By Proposition 
14.61 C is not /^-avoiding. 

Case 6. Assume (C \ v)* = {{}) . Then C\v = Am+n+i and so C = A{Dm,n) 

with V as the added vertex. □ 

Now we have all the necessary tools to prove Theorem 16.31 

Proof of Theorem \6.‘A By Proposition 13. 2Ul (IT!) implies ([2]). Each nuclear complex can be 
obtained by applying the unimodularity-preserving operations from Section [3]to A„U A^, 
which is unimodular by Proposition 13.211 Therefore ([3|) implies ([T]). We now show fl2l) 
implies ([3]). 

Let C be a /9-avoiding complex on q vertices. We show that C is nuclear by induction 
on q. For the base case, note that all simplicial complexes on 2 or fewer vertices are both 
/9-avoiding and nuclear. We may assume that C has no ghost vertices since adding and 
removing ghost vertices does not affect the properties of being nuclear nor /9-avoiding. In 
light of Propositions 14.41 and 15.41 and Lemma 13751 we only need to consider the cases where 
the 1-skeleton of C is or the union of two complete graphs, glued along a common 
















20 


DANIEL IRVING BERNSTEIN AND SETH SULLIVANT 


nonempty clique. Therefore, we can choose a vertex u of C that is in a facet with every 
other vertex. Since C is /9-avoiding, so is link.y(C). By induction, linkt,(C) is also nuclear. 
There are hve main cases. 

Case 1. This is the case where hnki,(C) = A(P) for a nuclear T). Let F denote the big 
facet of C and let u be the vertex of C that is not in F. Then C = A(C \ u), and C\u is 
nuclear by induction. 

Case 2. This is the case where hnk„(C) has ghost vertices. Since v connects to all 
other vertices, this is impossible. 

Case 3. This is the case where link.i;(C) = cone^(Am U A„) with p,m,n> 0. We split 
this into the sub cases, the hrst where m, n > 1 and the second where m = 0. 

Case 3.1 This is the sub case where m,n > 1. Proposition 13.91 gives C* \ v = 
cone^(Zlm,n) and so C* has Dm,n induced. Proposition 14.51 then implies that C* is ei¬ 
ther or coneP{Dm,n) with u as a ghost vertex. In the hrst case, C = 

cone^’''^(Am U A„). In the second case, C = A(cone^(Am U A„)). 

Case 3.2 This is the sub case where m = 0. Proposition 13.91 gives C*\v = cone^(Zlo,n) 
which has a ghost vertex u. Then C* \m is nuclear by induction and since u is still a ghost 
vertex in C*, this implies C* is nuclear. Proposition 16.21 implies that C is nuclear. 

Case 4. This is the case where link^(C) = cone^(Zlm,n) with p > 0 and m, n > 1. By 
induction, all proper minors of C are nuclear. Then Lemma 16.81 implies that C is nuclear. 

Case 5. This is the case where linki,(C) = A^ with k > —2. Since hnk^(C) has no 
ghost vertices k = q — 2. So C = Aq_i. □ 

7. Unimodularity of Non-Binary Hierarchical Models 

In this section we consider the unimodularity of hierarchical models that are non¬ 
binary; i.e. determining when .4^ d is unimodular for arbitrary d. We begin by showing 
how Theorem 16.31 makes this problem more tractable. Unlike in the binary case, the 
Alexander dual does not always preserve unimodularity. For example, .4.c,d is unimodular 
when C = AiUAi and d = (2, 2, 2, 3). However, .4c*,d is not (see Proposition l7.9p . Because 
of this, our proof of Proposition 13.101 is not valid in the non-binary case. The proposition 
statement does however generalize as we show with Corollary 17.41 Proposition 17.91 gives 
a list of known minimal pairs (C, d) with C nuclear where .4c,d is not unimodular. We 
then use this list to pose a simple question whose answer would be a step in the general 
classihcation. 

Proposition 7.1. Assume Ac,d D unimodular. Then for all d' such that d' < d compo¬ 
nentwise, Ac,A' is unimodular. 

Proof. Note that .4c,d' can be realized as a subset of columns of .4c,d- D 

The special case where d' = 2 gives us the following useful corollary. 

Corollary 7.2. If Ac ,a is unimodular, then C is nuclear. 

So in order to classify all unimodular hierarchical models we only need to consider nu¬ 
clear complexes. Therefore we can approach the general classihcation problem by looking 
at each nuclear C and identifying the minimal values of d that give rise to non-unimodular 
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Ac,d- Before proceeding with this, we show that the class of unimodular Ac,d is still closed 
under taking minors of C. However Corollary 13.71 fails in the non-binary case; we cannot 
freely take the Alexander dual of C. Because of this, our proof of Corollary 13.101 is not 
valid in the non-binary case. We will give an alternate proof that the class of unimodular 
^c.d is closed under taking links in C. We start with a useful proposition. 

Proposition 7.3. Let A G be a unimodular matrix with columns Assume 

is nonzero. Let A' be the matrix that results when we project A onto the hyperplane 
orthogonal to a„. Then A' is unimodular. 

Proof. We may assume A has rank r since otherwise we can delete unnecessary rows. Let 
B E have rank n — r such that = 0. Denote the columns of B as {bi}r=i- 

By Proposition 13.51 B is also unimodular. Denote the columns of A' as Then 



Let B' be the matrix with columns Then B' has rank n — roin — r — 1. The 

second case implies that b„ is a coloop in the matroid underlying B. Since the matroids 
underlying A and B are duals, this would imply that a„ = 0. But then A' = A which we 
assumed to be unimodular. So we can assume that rank(H') = rank(i?) = n — r and that 
7 ^ 0. Therefore rank(A') = r — 1 and so the dimension of its kernel is n — r, which is 
the rank of B'. We claim that = 0. From this it follows by Proposition 13.51 that 

A’ is unimodular. 

Now we prove the claim. Let b* be a column of B' (so 1 < i < n — 1). Letting hji 
denote the jth entry of bj, we have 


A'bf = 




Note that bjn = 0, so we can add the nth term to the above sum. This 

enables us to break it up as follows 


i=i 


E", 

i=i 


a.j. a, 




i)^n 


The hrst term is AhJ = 0. The second term is (^ ya^a a^Abf j ■ a„ = 0. So the claim is 
proven. □ 


Corollary 7.4. Assume Ac,d is unimodular. Let v be a vertex of C and let d' denote the 
vector obtained by deleting the entry for v from d. Then A \\r^v ^.c d' is unimodular. 
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Proof. Let A = ^iink„c,d' and let B = By the remark following Lemma 2.2 in Q 

and the lemma itself, we can write 


'A 

0 . 

. O' 

0 

A . 

. 0 

0 

0 . 

. A 

_B 

B . 

■ B_ 


Assnme A G and B G R^^"' and so ^ By Proposition I2.9fl3l) . we may 

remove rows from A to make A, B, A all have fnll row rank. So assnme that they do. Let 
A^ denote the matrix that resnlts when we project Ac,d onto the snbspace orthogonal to 
the last {d — l)n colnmns. Then if 1 < i < n, the ith colnmn of A! can be expressed as 
the ith colnmn of Ac,d minns a linear combination of the last {d — l)n columns of Ac,d^ 
all of which are 0 in the top m rows. So this means that the top m rows of A! are 


{A 0 ... 0) . 


Furthermore, since A and v4.c,d both have full row rank, the hnal {d — l)n columns of 
Ac,d have rank [d — l)m + 1. Therefore A! has rank m. Since A also has rank m, we may 
delete the bottom (d — l)m + I rows of A! without affecting the rowspace and therefore 
unimodularity. So the matrix (A 0 ... O) , and therefore A, is unimodular. □ 


Our proof of Proposition 13.11 is still valid in the non-binary case. Therefore, we have 
the following. 


Proposition 7.5. Assume C is a minor of C. If Ac,d is unimodular, then so is .4,c',d' 
where d' is the restriction of d to the vertices that are in C . 

The proofs of Propositions 13.121 13.171 generalize to the non-binary setting and so we 
have the following. 


Proposition 7.6. If the pair (C, d) gives rise to unimodular Ac,d, so do (cone^(C), d') 
and {G{C), d") where d' = (d ci ... Cp) and d" = (d c) for any c,Ci,...,Cp> 2. 


Proof. In the case of a cone, the matrix AconeP{c),d' is a block diagonal matrix: 




(•Ac,d 
0 


conef (C),d' 


0 

‘^C,d 


0 


Vo 0 ■■■ .Ac,d/ 


with C 1 C 2 ■ ■ ■ Cp blocks down the diagonal. In the case of adding a ghost vertex, the matrix 
•A.G{c),d' has repeated columns: 

•A.GiC),d" = (-^C.d •A.C,d ■ ■ ■ -^C.d) 

with c blocks. □ 


We can also extend Proposition 13.211 
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Proposition 7.7. Let C = AmU and d G Let M,N denote the vertex sets of 

C in the A^, A„ respectively. Let V be the complex with facets {1, 2} and let e = (ci, 62 ) 
where 

ei = n 62 = JJ 4. 

v&M vGN 

Then Ac,d = Av,e- This matrix is unimodular. 

Proof. We can see that Ac,d = by inspection. The matrix Av,(ei,e 2 ) ^be vertex 
edge incidence matrix of a complete bipartite graph with ci and 62 vertices in each set of 
the partition. Snch vertex-edge incidence matrices are examples of network matrices and 
are hence totally nnimodular 0 Ch. 19]. □ 

We also note that Proposition 13.151 holds in a slightly more general setting. 

Proposition 7.8. Let (C, d) he such that .4,c,d is unimodular. Then .4,A(c),d' is also 
unimodular if d' = (d 2 ), 

Proof. The proof here is similiar to the proof for Proposition 13.151 Note that A(.Ac,d) = 
^A(c),d'- Hence Theorem 7.1 in jl^ implies the proposition. □ 

Propositions 17. II and 17.51 suggest that we might be able to use a short list of pairs (C, d) 
that give rise to non-unimodular Ac,d to eliminate the remaining non-unimodular pairs. 
The following proposition lists all non-binary minimally non-unimodular complexes that 
we presently know. 

Proposition 7.9. The following pairs give rise to non-unimodular Ac,d 

(1) C = A(Ao U Aq) with facets {12, 23,13} and d = (3, 3, 3) 

( 2 ) C = Dl l with facets { 12 , 23, 34,14} and d = ( 2 , 2 , 2 , 3) 

(3) C = A(G(Ao U Ao)) with facets {12,13, 234} and d = (4, 2 , 2 , 2 ) 

(4) C = A(T)i,i) with facets {1234,125, 235, 345,145} and d = ( 2 , 2 , 2 , 2 , 3). 

Proof. By Table 1 in [^, ([T]) is not unimodular. Computing a Graver basis with 4ti2 [l[ 
shows that ([2]) and (jl]) are not unimodular. Computations are available on our website 
[^. If C and d are as in ([3]) then Ac d = Aj* 2 which is not unimodular by Proposition 

O ’1’ ^ 

From these non-unimodular examples and the constructive results earlier in the section, 
we can classify all unimodular hierarchical models that are nuclear with nucleus Dm,n, 
m,n > 1. 

Theorem 7.10. Let C be a nuclear complex with nucleus Dm,n with m,n > 1. Then Ac,d 
is unimodular if and only if 

( 1 ) = 2 for all v in the original Dm,n 

(2) = 2 for all v corresponding to a Lawrence lifting. 

Proof. The fact that all such Ac,d are unimodular follows from applying Propositions 17.61 
and 17.81 We now show that these are the only unimodular Ac^d with this type of complex 

C. 
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Note that since taking cone vertices commntes with adding ghost vertices and Lawrence 
liftings and this always preserves nnimodnlarity, we can assnme there are no cone vertices. 
Now if C is a nnclear complex obtained by snccessively adding ghost vertices and Lawrence 
vertices to then the induced subcomplex on just the Dm,n and Lawrence vertices is 

an iterated Lawrence lifting of D^ n- If any of the vertices of v that are Lawrence vertices 
have > 2, then by taking suitable links, one obtains the complex isomorphic to item 
(4) in Proposition 17.91 On the other hand, if > 2 for some v in the original Dm,ni then 
by taking suitable links, one obtains the complex isomorphic to item (2) in Proposition 
17.91 In summary, if the conditions of Theorem 17.101 are not satisfied, there is a minor of 
C that is not unimodular. □ 

It remains to consider the nnimodnlarity of ^c,d when C is a nuclear complex with 
nucleus AmU A„. We do not have a general characterization in this case. By Propositions 
ESI and 13.181 we may assume that C has no cone vertices. By Proposition 17.71 we may 
assume that m, n = 0. Note that A(Ao U Aq) = CIA 2 , and our results above (or [^) show 
that this is unimodular if and only if some dy = 2. Similarly, the iterated Lawrence liftings 
A^(Ao U Aq) = cIAp_|_i is unimodular if and only if at most two of the dy are greater than 
2 . 

We begin moving into uncharted territory when we consider cases where ghost vertices 
are involved. The simplest case to consider is the complex A(G'(Ao U Aq)) which is the 
complex on {1,2, 3,4} with facets {12,13,234}. Here 1 is the Lawrence vertex and 4 is 
the ghost vertex. By Propositions 17.51 and 17.91 1|3|), we must have di <3. If di = 2 then 
Proposition ^. 81 implies that v4.c,d will be unimodular. So we only need to consider the case 
where di = 3. Propositions 17.1117.51 and 17.91 ([T]) imply that one of ^2 or ds is 2, so assume 
d 2 = 2. This leaves us with the following question, which would classify nnimodnlarity 
forC = A(G'(AoUAo)). 

Question 7.11. Let C = A(G(Ao U Aq)) be the complex on {1,2,3,4} with facets 
{12,13,234}. Let d = (3,2,d3,d4). For what values of d 3 ,d 4 is the matrix Ac,d uni¬ 
modular? 
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